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Abstract 

In this paper we make further refinements to the duality proposed between N — 1 SQCD 
and certain string (supergravity plus branes) backgrounds, working in the regime of com- 
parable large number of colors and flavors. Using the string theory solutions, we predict 
different field theory observables and phenomena like Seiberg duality, gauge coupling and its 
running, the behavior of Wilson and 't Hooft loops, anomalous dimensions of the quark su- 
perfields, quartic superpotential coupling and its running, continuous and discrete anomaly 
matching. We also give evidence for the smooth interpolation between higgsed and confining 
vacua. We provide several matchings between field theory and string theory computations. 
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1 Introduction 

Maldacena proposed the AdS/CFT Conjecture [T| almost ten years ago and since then, 
great advances have taken place to refine and extend the original duality (see [2], [3] and 
references to those articles). One of the lines of research that captured the interest of many 
physicists is the use of the ideas and framework presented in [TJ-[3J, to formulate duals to 
phenomenologically relevant field theories. The final aim of this line is to find a dual to 
large N c QCD, that according to early ideas of t'Hooft [I] should be described by a theory 
of strings. 

A slightly less ambitious project (and perhaps a necessary step to overtake the case of 
QCD) is to find duals to phenomenologically viable field theories with minimal SUSY. Again, 
remarkable advances on this area of research have been achieved in the last nine years, see 
for example [5] and references to those works. One well-known limitation of the papers 
mentioned in [5] and successive work on that line is the fact that the dual field theory contains 
fields transforming only in the adjoint representation (also in bi-fundamental representations 
for the cascading theories), hence they could be thought of as SUSY versions of Yang- Mills. 
The rich dynamics added by the presence of fields transforming non-trivially under the center 
of the gauge group, was absent in the examples of [5]. 

To ameliorate this situation the introduction of 'flavor' branes (that is spacetime filling 
branes usually with higher dimensionality than the 'color' branes) was considered. Indeed, 
the flavor degrees of freedom need the presence of an open string sector, and this can only be 
achieved by placing in the background Nf flavor branes, on which the SU(Nf) symmetry will 
be realized. This problem was first treated in the 'probe' approximation, where the small 
number of Nf flavor branes was not changing the geometry generated by the large number 
of N c color branes. This non-backreacting approximation has a diagrammatic counterpart 
in the dual QFT, called the quenched approximation. The basic idea is that when funda- 
mentals run on internal loops of the diagram, the putative diagram is weighted by powers 
of Nf/N c and therefore suppressed in this approximation. In other words, the quenched 
approximation boils down to considering the path integral of a theory with adjoints and fun- 
damentals and approximating it by the path integral of the theory with adjoints only, where 
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the fundamentals are only external states (or run on external lines). The approximation is 
good if the quarks are very massive, hence difficult to pair produce, or if Nf/N c —>■ and 
this is the non-backreacting condition mentioned above. The papers that shaped this line of 
research are listed in [6] and there were interesting follow-up works. The outcome is a set of 
spectra of mesons for different QFT's duals, a nice picture of chiral symmetry breaking and 
phase transitions. 

A natural question then is if some relevant dynamics is lost when working in the t'Hooft 
expansion iV c — > oo, Nf = fixed. Indeed, one can see that this is the case. From experience 
in lattice field theory, the quenched approximation is good to compute static properties (like 
the spectrum of QCD), but works poorly when trying to address thermodynamical, phase 
transitions or finite density problems (see for example [7]). Obviously quenching will be a bad 
approximation if the quarks are light and/or if the number of fundamentals is comparable 
to the number of adjoints. 

A natural problem is then to understand how to backreact with the flavor branes or, in the 
dual field theory, how to 'un-quench' the effects of the fundamentals. From the diagrammatic 
viewpoint, this unquenching appears when considering the Veneziano Topological Expansion 
[8] instead of the t'Hooft expansion [3]. 

On the string side, this problem did not go unnoticed and the approach proposed was to 
find solutions to an action that puts together the supergravity action (for the color branes) 
and the Born-Infeld-Wess-Zumino action (for the flavor branes). 

S = S 'gravity + Sbranes = Sjjb/A + SbI + Sy/Z (1-1) 

Cases where the flavor branes dynamics, represented by Sbranes in eq. fll.ip are localized 
by some delta-function (that is, finding solutions in pure supergravity for intersecting D p - 
-D p+ 4 branes), have been discussed, see for example [9|. In these cases, typically there is 
a curvature singularity at the position of the delta-function, which may hide some of the 
interesting dynamics arising from Sbranes- 

On the other hand, cases where Sbranes has support in the whole spacetime were first 
studied in non-critical string set-ups [10] and later in type II supergravity where a variety 
of examples were constructed, ranging from duals to a version of minimally Supersymmetric 
QCD (SQCD) pj] and [32], N = 2 version of SQCD [13] and the addition of fundamentals to 
cascading field theories/backgrounds [TJ]. Recently, a localized solution for the D4 — D8 — DS 
system was found in [15] . Some elaborations on the thermodynamics of the N = 1 SQCD 
example mentioned above can be found in [16] and applications to SUSY breaking in [17]. 
From a more formal point of view it was shown in [18] that our intuitive way of adding flavor 
branes can be put on a formal basis where SUSY is guaranteed whenever the flavor branes 
are calibrated, and the second-order Einstein equations are automatically satisfied once the 
first-order BPS and Maxwell equations and Bianchi identities with sources are solved. 

Notice that in the context of critical string theory, backreacting with the flavor branes as 
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indicated in eq. fll.ip is a natural continuation of the approach in [6] and indeed the correct 
way to introduce the fundamental degrees of freedom with SU(Nj) global symmetry. 

In this paper, we will study the dual to minimally SUSY QCD that was proposed by the 
present authors in [TT] . Indeed, in [TT] a set of solutions to type IIB supergravity plus branes 
has been found and a large number of checks presented, where it was shown that the type 
IIB backgrounds were capturing non-perturbative effects of this version of SQCD. 

We will present a new set of solutions to type IIB supergravity plus branes, dual to the 
field theory in different vacua and will provide a large number of matchings to show that 
our new solutions indeed capture new qualitative and quantitative information of the dual 
QFT. 

The organization of this paper is as follows: in section [2] we will review some material 
from [IT] , set up the notation for the string backgrounds and dual SQCD. We also discuss 
some aspects of the SQCD-like theory that can be learned from the background without 
knowing the explicit solution, like Seiberg duality and the matching of R-symmetry anoma- 
lies, and we discuss some aspects of the field theory that will be relevant for the rest of 
the paper. In section [3] we present the new solutions that motivate this paper. We show 
that an expansion for large values of the radial coordinate (the UV of the dual SQCD) 
can be smoothly connected with the IR expansion of the solution. We discuss in detail 
the numerics and present plots to clarify ideas. In section @] we study the predictions that 
these new solutions make for the physics of the SQCD theory; they include definitions of 
gauge coupling, anomalous dimensions, quartic coupling and computation of the beta func- 
tion, continuous and discrete anomaly matching, potential between quark-antiquark and 
monopole-antimonopole pairs as computed by Wilson and 't Hooft loops and different as- 
pects of the QFT that we learn from the string solution. We close with some conclusions 
and future work to appear. We complement our presentation with an appendix on a nice 
technical subtlety. 

2 Brief Review of the dual to SQCD 

In this section, we will go over some points of [TT] that will be relevant in the rest of the 
paper. First, we will review the two types of backgrounds that follow from the work 
After that, we will discuss some aspects of the dual QFT. 

2.1 The two types of backgrounds 

We presented in [TT] a string dual to a version of iV = 1 SQCD. The construction proceeded 
as explained around eq. fll.ip by adding flavor branes to a given supergravity solution and 
finding a new background where the effects of flavor branes were encoded. The unflavored 
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solution we based our construction on, was proposed to be the dual to N = 1 SYM in [19] 
(the solution was originally found in 4d gauged supergravity in [20]). We considered the 
backreaction of the flavor branes, that in this particular case are Nf D5 branes wrapping 
a non-compact two-cycle. We found two types of solutions, which in the following we will 
refer to as solutions of type A and solutions of type N. 

Backgrounds of type A look quite simple, consist on a metric, a Ramond three-form, a 
dilaton and they read, in Einstein frame and setting a' — g s — 1 (we have slightly changed 
the notation with respect to [TT]). 



ds 2 



(3) 



■Hp) 

e 2 



dxl 3 + AY(p)dp 2 + H(p)(d9 2 + sin 2 9dip 2 ) + G(p)(d6 2 + sin 2 9dip 2 ) 
+Y(p) (dip + cos dd(p + cos 9dip) 2 ^ 



— sin 6d0 A dip + 
A 



N f -N c 



sin OdO A dip A (dip + cos 9dip + cos 6 dip), 



(2.1) 



where we have used the coordinates x = (x^, p, 8, ip, 8, ip, ip). The functions in eq. ()2.ip . satisfy 
a set of BPS eqs (that solve the Einstein, Maxwell and Bianchi eqs). The derivation of these 
BPS eqs was discussed in detail in [llj. They read: 



H' = -(N c -N f ) + 2Y 
G' = -— + 2Y 

2 V f H 2 G 

(N e - N f ) iYc 
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2Y 2 
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I, 
G' 



AY 



(2.2) 
(2.3) 
(2.4) 
(2.5) 



On the other hand, solutions of type N are slightly more involved and can be seen as a 
generalization of the type A backgrounds. The dilaton is a function of the radial coordinate 
<p = <f)(p), while the metric and Ramond three form read, 



ds 2 



tip) 
e 2 

p2<?0) 



+ 



dx* t3 + e 2k(p) dp 2 + e 2h(p) (d8 2 + sin 2 ddip 2 ) + 
((u>i + a(p)dO) 2 + (uj 2 - a(p) sinOdip) 2 ) + 



(u 3 + cos 6 dip)' 



N c 

F (3) = T 



(Cji + b(p)d6) A (uo 2 - b(p) sin 9 dip) A (uj 3 + cos 9 dip) + 



b'dp A (—d9 A lui + sin 9 dip A Cu 2 ) + (1 - KpY) sin 9d9 A dip A a; 3 
sin 9d9 A dip A (dip + cos 9dip) . 



(2.6) 
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Where Ui are the left-invariant forms of 577(2) 

Ox = cos ipdO + sin ip sin 6d(p , 
cD 2 = — sin ipd9 + cos ^ sin 6d(p , 

ujz = dip + cos8d(f . (2-7) 

For the case of backgrounds of type N we also found a set of BPS eqs and solutions to it, 
by expanding in series near p = (the IR of the dual QFT) and p — > 00 (the UV of the dual 
QFT) and then showing that the two expansions can be numerically matched in a smooth 
way. It was for the backgrounds of type N that we presented a considerable number of 
checks showing how they captured non-perturbative aspects of SQCD [TT] . 

In this paper, we will focus on backgrounds of type A above, finding a new set of solutions 
for the functions H, G, Y, <p and we will study the strong-coupling effects that they predict 
for the dual QFT. This kind of solutions was first overlooked in [TT] , because at that time 
we could only find badly singular type A solutions. We believe the new set of backgrounds 
we derive in the present paper, which are characterized by a good IR singularity, nicely 
complement the overall picture. 

After having set up the stage, it is now convenient to discuss some subtleties about the 
dual theory, that as we anticipated above, is a version of N = 1 SQCD. 

2.2 The proposed dual field theory 

As it is well known, all duals to non-conformal gauge field theories relying on the supergravity 
approximation like the ones in [5] , [19] or any other model available, are characterized by the 
presence of extra modes or some other kind of UV completion 0. This UV completion is not 
necessarily an ugly feature, one could perhaps be interested in the field theory with the UV 
completion that string theory provides (like in the case of cascading theories) or perhaps, 
there are ways to distinguish when a given field theory correlator is affected or not by the 
extra modes. In the particular case of interest in this paper, the field theory dual to the 
background in [19] [20J was proposed in the paper [19] to be N = 1 SYM plus an infinite 
set of KK modes, which is just another way of saying that the UV completion is a higher 
dimensional theory. The dynamics of these extra modes was studied in detail in [21] and a 
way to disentangle the effects of the UV completion from the Super- Yang-Mills modes was 
proposed in [22]. In summary, the picture of the dual field theory of [19] is clear by now. 
Simplifying things a bit, one could say that the field theory without flavors is described by 

lr This can be avoided in the non-critical string approach, but in that case the supergravity approximation 
is not reliable. 
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a lagrangian of the form (in fields and superfields notation) 

L = -Tr[^F%, + iXDX - J2( D M 2 ~ M 2 KK <& + $ h {iD - M KK )^ k + V{<f> k , A)] = 

k 

= [ d4d J2 $* eV$ * + / d2 ° (w a W a + J2 + h.c. (2.8) 

k k 

where, schematically, we wrote a lagrangian describing N = 1 SYM coupled to an infinite 
tower of SUSY KK modes, with a given superpotential W($fc) (for a detailed expression 
see [21]). 

We proposed in [11] that the addition of the flavor branes introduces two chiral superfields 
Q, Q transforming in the fundamental and anti-fundamental of the SU (N c ) gauge group, 
respectively. The dynamics and the coupling of the quark superfields with the already 
existing adjoints was proposed to be of the form, 

L = J d A 6(Q ] e v Q + Q j e v Q) + J d 2 6Q$ k Q + h.c. (2.9) 

where is a generic adjoint field. So, putting together (12.81) and (12. 91) . we have a lagrangian 
that in superspace will schematically read, 

L= j d 4 e(Qie v Q + &e v Q + J2®l eV ®k) + I d 2 e(w a W a + ^Q$ fc Q + W($ fc )) + h.c. 

(2.10) 

Let us consider the simplest case in which the superpotential contains only a mass term 
W($fc) = Hk&k- The occasional presence of higher powers of the chiral multiplet of KK 
modes will be irrelevant for the IR dynamics as we will see below. Since we are interested in 
studying the field theory at low energies, we can integrate out the scalar KK modes. This 
should be taken with some care because the mass scale of the KK modes is comparable 
with the scale of strong coupling which we are interested in, so the integration-out of these 
additional modes is not totally clean. After integrating over the dynamics of the scalar 
multiplet, we get a dynamics for the quarks and vector superfields of the form 

L = Jd 4 6 (Q^c v Q + Q ] c v Q) + J d 2 9 (w a W a + i [Tr(QQ) 2 - ^(TrQQ) 2 ^ + 0(QQ 3 )) 

M C (2-11) 

So, we arrive to a theory with the field content of iV = 1 SQCD, but with an important 
dynamical difference: the flavor symmetry is SU(Nf). The root of this difference is basically 
that the 'mother theory' (12. lOf) has already SU(Nf) and not SU(Nf)i x SU(Nf) r because 
of the term Q^ k Q- This same qualitative difference from SQCD or QCD will generically 
appear when adding flavors to field theories that could be thought of as coming from an 
N = 2 theory. On the other hand, when color and flavor branes are codimension six, as for 
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instance in the last paper of ref. [6] and [23j[T5], this kind of coupling is not present. Another 
interesting point is that if the adjoint superfield transforms in U(N C ) instead of SU(N C ), 
then the double trace term in the superpotential is absent. 

In the following, only the quartic coupling for the fundamentals will be considered. In 
section 12.2.21 we will argue that higher order couplings are irrelevant. We present here a 
heuristic complementary argument which favors the appearance of the quartic coupling and 
not the higher ones upon integrating-out the adjoint scalars. The superpotential for the 
adjoint scalars is of the form [21], 

W($ fc ) =Tr[$ 1 [$ 2 ,$3] + ■ ( 2 - 12 ) 

i 

The "democratic" way in which the Nf flavor branes are treated suggests that the QQ matrix 
is proportional to the identity (notice that QQ here is not a meson operator, gauge indices 
are not contracted). When we integrate out the adjoints, and insert back their value in 
W($fc), the higher order contributions involving QQ vanish using conmutation relations. 

We would like to mention at this point that, in order to get our BPS eqs and their solu- 
tions, we introduced in [11] a smearing procedure. This has some effect on the superpotential, 
as discussed in that work. We will not insist here on the fact that qualitative and quanti- 
tative aspects of the localized solution are well captured by the smeared solution. For more 
discussions, we refer the interested reader to our forthcoming paper [24|. 

In the paper [TT] , backgrounds of type N in eq. CI2 . €31) above were shown to capture a variety 
of non-perturbative aspects for a field theory like (12.111) . In the rest of this paper, we will 
show how backgrounds of type A also encode interesting dynamics of ( 12. lip . We now turn 
to study some aspects of the theory defined in eq. (12.111) that will be relevant for the rest of 
the paper. 

2.2.1 A couple of important subtleties 

We want to explain briefly some nice subtleties that arise when considering the theory (12. lip . 
Let us start with the superpotential in (12.111) . 



W tree = - 
[I 



Tr(QQ 2 ) - ^-{TrQQf (2.13) 



N, 



c 



It is known that when analyzing the space of vacua for this theory the quantum mechanical 
superpotential differs from the classical one in eq. (I2.13p . Indeed, if Nf < N c an Affleck- 
Dine-Seiberg superpotential [25] is induced 

w = w '- + ^- N 'i^m) ' (2 - 14) 
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which changes qualitatively the dynamics, for example, when looking for solutions to the 
F-term eqs, not allowing the (diagonalized) meson field M = QQ to have zero eigenvalues. 

When Nf > N c a different quantum superpotential is induced. Actually, for the cases 
Nf = N c , Nf = N c + 1 and Nf > N c , three different superpotentials are induced that present 
different space of vacua. Here, we will content ourselves with stating that for the case of 
Nf > N c , the field theory presents two possible behaviors, one in which the meson matrix 
has non-zero eigenvalues and another branch where the meson matrix has zero determinant 
(the meson matrix can always be diagonalized). For details, see (26], [27]. 

We will explain in section [3] how these two possible behaviors are matched by the super- 
gravity solutions of type A and type N. 

The second interesting subtlety refers to the large p behavior of the solutions we found. 
As we will see in section I3.4[ the careful study of the UV (large p expansion) of our solutions 
suggests that in the dual field theory an operator of dimension six is getting a VEV. Probably 
as a consequence of this, some well-known or expected behaviour of the IR dynamics of the 
theory (12. lip , see [29], will not be realized by our string background, that in turn is dual to 
the field theory but in a particular vacuum. Indeed, in the following we will check that the 
string solutions seem to agree with the interpretation of SQCD with quartic superpotential 
in the UV (large p) region but we will observe that in general the expected IR fixed point 
is not attained. Of course, this depends on the precise vacuum on which the theory is 
realized. As a matter of fact, if the theory flows to the IR Seiberg's fixed point, i.e., the 
quartic coupling flowed to zero, there is an enhancement of the flavor symmetry at the IR 
fixed point. It is hard to see how a brane construction of the kind we are using can reflect 
this fact. Thus, there must be something in the field theory dual to the presented solutions 
preventing the vanishing of this coupling. It would be very interesting to understand this 
point better. In section 14.61 we will comment on the actual IR behaviour which can be read 
from the geometry. 

2.2.2 Beta functions 

We know that in an N = 1 field theory, there are lots of simplifications for the beta functions 
of different couplings. For the gauge coupling, the well known NSVZ- Jones beta function 
rGcicls 

ft =-(£^)K -"/<!-•»>) < 2 - i5 > 

where 7q is the anomalous dimension of the quark superfields. In the following, we will 
be interested in the wilsonian beta function, for which the denominator in (I2.15P is absent. 
Suppose that the superpotential for our model is (here we neglect the double trace terms of 
the form -^-(TrQQ) 2 , that will not change the beta functions) 

W = k(QQ) 2 + A(QQ) 3 + z(QQ) 4 + ... (2.16) 
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Let us remind the standard result in four-dimensional SUSY gauge theories; for a superpo- 
tential of the form W = hQ n A m if the fields have (classical) dimensions [A] = mass a , [Q] = 
mass q we can form the quantity = 3 — nq — ma and the adimensional coupling h = hfi~ Ah 
( where [i is some mass scale). Then, the beta function for the coupling h is 

~ ti th 

h = h[-A h + - lQ + - lA ] (2.17) 

where 74, 7q are the anomalous dimensions of each field. Let us apply this to our case ( 12. 16ft . 

A* = «[1 + ^7q], ^ = A[3 + | To ], p s = ~z[5 + ^ 7 q] (2-18) 
We will argue in section 14.1.11 that in the solutions for Nf < 2N C , near p — > 00, one finds: 

T0 = -i-e, (iV / <2iV c ) (2.19) 

with e some small correction in the UV (but at lower energies it may become large). Thus, 
when Nf < 2N C the quartic coupling k is relevant (becomes marginal at Nf = 2N C ) , but the 
others are irrelevant. When Nf > 2N C , we have 7q > — ~ and all couplings are irrelevant 
(their beta functions are all positive). So we see that dropping (QQ) 3 and higher orders 
from the lagrangian (12.111) is well motivated if we are interested in the IR dynamics. Also, 
the argument around eq.f l2.12p suggests that the higher order operators do not occur at all, 
so in the following we will consider the case of quartic superpotential only. In section I4.1.H 
we will show how solutions of type A reproduce the beta functions for the gauge coupling 
in (12.15j) and the anomalous dimensions in ( 12. 19ft . 

Let us now focus on another interesting property of the field theory ( 12. lip . Seiberg duality. 

2.3 Seiberg Duality 

A well known property of N = 1 SQCD-like theories is Seiberg duality [28]. In the particular 
case of field theories like the one in eq. (12. lip , this equivalence is specially nice as explained 
for example in [29J, where it was shown that the duality can be exact and not just an IR 
equivalence of two theories. In [TI], we explained how the duality proposed by Seiberg is 
manifest in our backgrounds of type N. Here, we will concentrate on showing how type A 
backgrounds reflect this duality. 

Indeed, without knowing the solution to the BPS system of equations (I2.2l) - ()2.5p . we can 
observe that the equations and hence their solutions are invariant under: 

H «-> G, A c -> N f -N c (2.20) 

with 0, Y, Nf unchanged. This change has the correct numerology to act as Seiberg duality. 
Indeed, we can see that the transformation (12.2011 is idempotent. 
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Let us comment on the correspondence between Seiberg duality and the transformation 
(I2.20p . The reader may want to think about it in this way: if we take backgrounds of the 
form (12.11) . characterized by BPS eqs. fl2.2p - fl2.5l) . then it is observed that a solution for the 
functions H, G, Y, <p, dual to a field theory with Nf flavors and N c colors, automatically 
produces 'another' solution for the functions G, H, Y, that we can associate with the field 
theory with Nf flavors and Nf — N c colors. The fact that the two solutions are the same 
solution (up to a coordinate redefinition), indicates that the two field theories are actually 
the same theory. This nicely matches with the exactness of Seiberg duality for theories like 
(12. lip , as discussed in detail in [29] . 

One natural question that may appear at this point is whether the Seiberg duality change 
(I2.20p is only a property of the SUSY system or, on the contrary, all solutions (SUSY and 
non-SUSY) to the eqs of motion that dictate the dynamics of type A backgrounds have the 
property above. To sort this out, notice that not only the BPS equations, but also the type 
A ansatz (12.11) is invariant under (I2.20p if we relabel: 

9^9, (p^fi (2.21) 

We can see that Seiberg duality is basically the interchange of the two-manifolds labelled 
by {9, ip) and (#, y5). Thus, even non-supersymmetric solutions of the resulting equations of 
motion, would have a Seiberg-dual counterpart. 

This kind of reasoning should nevertheless be taken carefully. Indeed, Seiberg proved [28] 
that in dual pairs there was matching of global anomalies and that deformations of moduli 
spaces could be put in correspondence. It is therefore necessary to study how anomalies for 
the R-symmetry are manifest in each theory in the dual pair. The prescription to compute 
R-symmetry anomalies was studied in [11] and, in following sections, we just apply that 
information for the matching of R-symmetry anomalies on both dual pairs. 

Summarizing, we propose the following interpretation: since a gravity solution corresponds 
to a field theory in a given vacuum, the fact that two different field theories correspond to the 
same gravity solution means that they are the same theory, they contain the same information 
(they are Seiberg dual). However, the way of extracting gauge theory information such as 
couplings or rank of the gauge group from the geometry can depend on the interpretation, 
on how one wants to understand -in view of the change (12.211) - the role of the different cycles 
in the geometry among other things. 

2.3.1 R-symmetry anomalies 

As usual for anomalies, we look into the Ramond fields, in this case C<i- One might be 
confused, since we write a potential C<i when eq. (12.11) implies that dF 3 ^ (that is the 
Bianchi identity indicates the presence of sources). The point that we made in [11] was that 
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one should first restrict the background to the cycle, 

'0 = 0, p = 2rr-^, ip), (2.22) 



since this is where color branes are wrapped around. On this sub-space F 3 is exact and one 
can write a potential C 2 . For the theory with Nf < 2N C , that we call e-theory (for electric), 
and for the Seiberg dual theory, with Nf > 2n c , that we call the m-theory (for magneticjj, 
the potentials are, 

N 

C 2e = (if>- ^0)^(2 - x)sin6d6 A dip 

C 2m = (V> -il) G )^{x-2)sinddd Ad<p (2.23) 

We have defined x = Nf /N c and x = Nf/n c as in [11]. Notice that both expressions in (12.231) 
are the same expressed in electric and magnetic variables respectively. 

Imposing that the instanton that is represented by a Euclidean Dl-brane wrapping the 
contractible two-cycle (I2.22|) is not observed in Bohm-Aharonov experiment [30], leads to 
some particular values of the ip angle (or selects some possible translations in ip, if we make 
ip — > ip + 2e). Indeed, remember that ip — > ip + 2e (with e in [0, 2ir]) was the change of 
R-symmetry in the unflavored case [19] and the unobservability of the shrinking instanton 
implies for the partition function of the putative Euclidean brane, 

Z = Z[0}e^f C2 -> e ieNc{2 - x) e^f C2 x Z[0] = Z (2.24) 

where Z[0] is the partition function for this probe including the Born-Infeld part; this in 
turn leads to, 

_ 2jm _ 2nk 

6e _ n AT AT ' e ' m ~ AT O \£.AO) 

2N C — Nf Nf — 2n c 

The relation in eq. (12.251) should be interpreted as the matching of the anomaly in the correla- 
tor of two gauge currents and the R-symmetry current, < J(R)J(N C )J(N C ) > on the theories 
of the dual pair. We can see that e e = e m as k = (1, ....Nf — 2n c ) and n = (1, ...2N C — Nf). 
This will work in the same way for the type N backgrounds, since it is the same RR field 
C2 that is relevant. This is not one of the diagrams alluded by the 't Hooft criteria, those 
will be discussed in section 14.81 

The fact that the R-symmetry is broken by anomalies to the discrete subgroup Z 2 n c -n s 
suggest that the R-charge for the quark superfield is R[Q] = |. Indeed, as we will see in 
section 14.1.11 and in section I4~8l the result of eq. (12.250 implies that the R-charge of the quark 
superfield is R[Q] = 1/2 (in the UV). The anomalous dimension of the quark superfields 
must then be (in the case of an UV fixed point) 7q = — |, in agreement with what we 
observed around eq. (12. 19ft . 



2 The relation is the usual n c = Nf — N c 
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2.4 A jewel of the 1970's 



Before leaving this field theory section, we would like to remind the readers about nice work 
done three decades ago. 

Several different suggestions have been made to address the problem of confinement. The 
idea of 't Hooft and Mandelstam [31 J was to think about confinement as a dual Meissner 
effect, where a 'vortex' of chromoelectric flux joined together two quarks giving an area law 
for the Wilson loop|§ The 'abelian-projection monopoles' proposed to be dual to the quarks 
[32] paved the way to interpret the phenomena of confinement in QCD as a dual Meissner 
effect. The following possible phases manifest in a non-abelian gauge theory: 

• Free phase: in analogy with QED; for example, if the number of flavors is very large 
in QCD. 

• Coulomb phase: the IR theory is conformal. 

• Higgs phase: here the monopoles are confined and quarks are free. 

• Wilson/Confining phase: here the quarks are confined and monopoles are free. 

• Oblique confined phase: quarks and monopoles are confined, dyons condense. 

In QCD, SQCD or in any QCD-like theory with fundamental matter, one can move smoothly 
(without phase transitions) between the different phases mentioned above [33]. We will see 
how the solutions of type A that we will find in section [3] realize these phases and their 
smooth interpolation. 

Let us now present the new solutions for the backgrounds of type A in eqs.f l2.1l) -f l2.5p . 

3 New solutions 

In this section we present the new solutions for backgrounds of type A described in eq. fl2.ip . 

Let us first explain why we are interested in solutions of type A. It is well known that even 
when singular, certain supergravity solutions can faithfully describe the dynamics of a dual 
gauge theory at strong coupling. The question of which singularities are allowed inspired a 
variety of criteria to discard those singular spacetimes that cannot be proposed to be dual 
to a given gauge theory and accept those that can. Among these criteria, we will choose 
the one in [31], that basically imposes on a singular spacetime to have finite value of gu at 
the singularity (in this discussion we are alluding to IR singularities, the UV singularities of 
these kind of theories are resolved in a stringy way). 

3 Needless to say, this regards the problem of confinement on Yang-Mills; in QCD, the presence of funda- 
mental matter allows to interpolate between the Higgs and the confining phase. 
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Different studies with supergravity solutions dual to field theories without flavors, have 
exemplified the correctness of this criterion. It is common lore by now (in the case of 
'unflavored' backgrounds), that solutions of the type A are singular with a so called "bad 
singularity" , that is, the putative background cannot be interpreted as being dual to a field 
theory, see for example [T9] . 

In our case, backgrounds of type A will have associated a value of gu = (in Einstein 
frame). If we analyze the BPS eq. for the dilaton- see the last eq of (12.51) - we immediately 
see that the dilaton has positive derivative if Nf > N c . Hence, type A backgrounds present 
a "good singularity" when the number of flavors is bigger or equal to the number of colors. 
This restriction very nicely matches what we explained in section 12.2.11 regarding the two 
possible vacua for Nf > N c . In other words @, backgrounds of type A describe N=l SQCD 
with the lagrangian (12.111) when the gaugino condensate and the VEVs of the meson matrix 
vanish. Conversely, we propose that backgrounds of type N studied in detail in [11] describe 
the same field theory, in a vacuum with non-vanishing meson VEV and gaugino condensate 
(we remind the reader that, in the unflavored case, it is well established that the function 
a(p) in (I2.6p is, roughly, the dual of the gaugino condensate and that type A solutions are 
obtained from type N by seeting a(p) = 0). 

In summary, for Nf > N c , type A backgrounds can be thought of as good duals to our 
N — 1 SQCD of eq. (|2.1ip in a vacuum with vanishing meson VEV and gaugino condensate. 
To describe these new solutions, we now present in detail the integration of the BPS eqs. (12.21) - 
(12. 5p for the type A backgrounds (12. ip . 

We start by noting that the difference H — G can be integrated immediately, 



G = H + 



N f - 2N C 



p-C 



(3.1) 



where C is a constant of integration. The function Y can also be written in terms of H, 



y = \d P H+\{N f -N L ) 



(3.2) 



This allows us to write the system as a single second order differential equation for H, 



d 2 H 
dp 2 



\d P H+\{N f -N c 



_ 2 dpH + N, 



H 



N f + 2d p H 
H+^^p-C 



(3.3) 



and the dilaton that can be obtained as, 



<t>{p) = 0o + 



N f -Nr. N r i 



AH 



+ AG 



dp 



(3.4) 



There is of course, one simple solution for Nf = 2N C (the case Nf = 2N C is the "conformal 
solution" studied in pi]) that reads H = N c /£, G = N C /(A - £), Y = N c /A, but even in 



J This interpretation was developed in discussions with Francesco Bigazzi and Aldo Cotrone. 
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the particular case Nf = 2N C there are more solutions that we will not discuss here [35J. 
Unfortunately we did not succeed in finding other exact solutions to (13.31) . Nonetheless we 
derive below the asymptotic expansions, near the UV and the IR, of the solutions and then 
show that a smooth interpolation between them exists. 



3.1 UV expansion (p — ► oc) 

As explained above, we are interested in solutions to (I2.2l) - (l2.5p with asymptotically linear 
dilaton. There are two distinct possibilities, depending on whether Nf is bigger or smaller 
than 2N C . 



3.1.1 N f < 2N C 

There is a one-parameter family (depending on the constant C that enters (13.11) ) of UV 
(large p) asymptotic solutions. For the different functions we have: 

Nf\ / N c \ N c Nf 1 

N, N c Nf N c Nf{N f -2N c -AC) 2 

A lQ(2N c -N f ) P 32(2iV c — Nf) 2 P 

N c N c Nf N c N f {Nf-2N c -AC) 3 

A OO/OAf AT" \ ' OO/OAf A 7" \ 2 ' * * ' 



4 32(2A^ C - Nf) r 32(2iV c - N fi 

_ x 8C + 2N c + 2N f _ 2 



d >* = 1 - 4 "" + 16(2^-^)^ + - (3 - 5) 

This coincides with the asymptotic expansion found in equation (4.19) of [H] for the solutions 
of type N, except for the fact that here there is an extra free parameter C, which in those 
solutions seems to be fixed to C = (Nf — 2N C ) / A. 



3.1.2 N f > 2N C 

Here again a one-parameter family of solutions depending on C is found, 

N f -N c Nf(N f -N c ) Nf(Nf-N c )(AC + 2N c -N f ) 

A lQ(N f -2N c ) P 32(N f -2N c ) 2 P 

N f -2N C N f -N c -AC Nf(Nf — N c ) _ x 

G = — pH - 1 — -p +... 

2 H A 16(N f -2N c y 

Nf - N c Nf(N f - N e ) 2 Nf(N f - N C )(AC + 2N C - N f ) , 

4 32(N f -2N c y 32(N f -2N c ) 2 9 

o , i -8C - 2V C + 3JV/ 2 . , 

= '- 4 " + 16 W -iWJ " + - (3 ' 6) 
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which coincides with the expansion for the type N backgrounds found in eq.(4.20) of [TT] if 
we set C = (Nf — 2iV c )/4. 



3.2 IR behaviour 

Let us now present the solutions to the BPS eqs. at small p, as they appear from solving 
(13. 3p together with (13. ip . (I3.2p and (13. 4p . We remind the reader that, as discussed in the 
beginning of section El we require Nf > N c so it is guaranteed that the quantity g u ~ e 2 
does not diverge in the IR. The space is singular at the origin, but we will call it a good 
singularity and proceed assuming that correct physics predictions for the IR of our theory 
of interest (12. lip can be obtained by doing computations with the type A backgrounds. 

There is quite a rich structure since we actually found three possible IR expansions that 
can be smoothly connected with the UV expansions of section 13.11 The three possible IR 
behaviors will be labelled as type I, type II and type III. Let us study them in turn. 



3.2.1 Type I expansion 

An obvious solution of (13.31) is H = \{Nf — iV c )(ci — p), where C\ is an integration constant. 
But eq. (13. 2p would yield Y = 0. Thus, this simple solution cannot be physical. Nevertheless, 
there are solutions that asymptotically approach this behaviour in the IR (as p — > — oo). 
The deviation from the simple linear solution is given by exponentially suppressed terms (as 
p — > — oo) in the following way, 

H(p) = ^y^ci ~P)+Y1 ^k(p)e 4kp (p - -oo) , (3.7) 

k>\ 



The Vk{p) are order k + 1 polynomials in p. The first one is 



V x = c 3 



P 2 ~ P(ci + c 2 + -) + CiC 2 + -(ex + c 2 ) + - 



(3.8) 



and the higher ones can be obtained iteratively. Notice that one of the integration constants 
can be reabsorbed by shifting the origin of p and that there is a relation among the different 
constants C = y(A^ — N c ) — yiV c . The leading IR behaviour of the rest of the functions 
can be readily obtained 

G(p) = ^{c 2 -p) + 0{e^) 

Y{p) = 4c 3 e 4p ( Cl -p)(c 2 -p) + 0(e 8p ) {p - -oo) 



e 



00 



, + O(e^) (3.9) 

V(ci - p)(c 2 - p) 

This IR asymptotics can be numerically connected to the UV behaviours described in section 
(13. ip . We depict some sample plots in section l3~3l 
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3.2.2 Type II expansions 



Differently from the type I expansion we presented in the previous subsection, equation (13.31) 
also admits solutions where H, G, or both vanish at a finite value of p, which we take to 
be p = without loss of generality. The way H and G reach the origin depends essentially 
on the value of the constant C in ( 13.31) . When C is strictly negative H goes to zero like pa 
whereas G goes to a constant. When C is strictly positive, H and G exchange roles: H goes 
to a constant while the leading behavior of G is pa. We call both these kinds of solutions 
type II. When C — 0, instead, H and G vanish at the same time, and their leading term goes 
like ps; we call this solution type III and postpone its presentation until the next subsection. 
More explicitly, for C < the following expansions solve equations (13.11) - (13.41) 

H =l'U>-+[^+N c -N f )p + h l c - ^pa+... 

hj N f \ 72C 2 + 20h 2 + 3C(10iV c - 7N 



G = -C + hip* + - -^J p + /ii ^> + . . . 

^$4(f + ^-^)^^ 72C2+20ft '^ (10JV '- 7JVj) ^- 



Nf — N c i 3G(N / -iV c ) 2 -/i 2 (2iV c + ;V 



whereas for C > 0, we have the following power series 

rr _ , i (h\ N f \ 72C 2 + 20h 2 + 3G(10iV c - 3iV» s 

, , /i? Ar \ , 72C 2 + 20/i 2 + 3C(10iV c - 3iV/) a 
G = h lP *- [-L + NAp + h c - ^pa + 



/ii l//i 2 , N c \ , u 72C 2 + 20h 2 + 3C(lON c -3N f) 



Y = ^-~2{3d + T) +k i ~^G^ ^ + "" 

, , N c i , 3C + hl(3N f -2N c ) 

= 0o+ 2/^ + Y2Ch\ P + - (3 - U) 

The constants hi, C and 0o are arbitrary. In particular 0o can be fixed to any value. For 
what concerns C and h%, instead, what happens is more interesting. While any choice of 
them gives a solution to equations (13. ip - (13 ,4p . requiring that a particular solution behaves 
in the UV as the functions in section 13.11 imposes an additional condition that enforces a 
relation between the useful values of C and h\. Figured] is the plot of such a relation for 
two fixed values of jf. All other possible choices still give a solution to ( 13. ip - ( 13. 4(1 . but the 
functions either reach a singularity at some finite value of p or reach the UV with a different 
asymptotical behaviour, in which we are not interested here (analogous to the one studied 
in section 8 of [TT]). 
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hi hi 




Figure 1: We represent here the relationship between C and hi that ensures that the evolution 
of a given type II expansion goes asymptotically to the behaviour we presented in section 
13.11 On the left we are considering negative values of C, corresponding to expansion (13.101) 
while on the right positive values are considered, which correspond to (13.111) . Dashed lines 
correspond to -J- = 1.5 while continuous lines correspond to = 2.5. 



3.2.3 Type III expansion 

When the constant C in (13.11) vanishes, type II expansions are not allowed anymore, but a 
new branch of solutions opens up, which we will call type III. In this case both H and G 
vanish at the origin of space. The expansion reads as follows 

i 1 2 
H = h lP * + j^{5N c - 7N f )p + -hip* + ... 

G = h lP $ - ^(5iV c + 2N f )p + h lP l + ... 

h x _ 2 N f 4, i 

Y = — p 3 l 4- -hxp-i + . . . 

6 H 10 9 H 

3Nf 2 10N*-10N c N f + 7N 2 f 4 

* = *o+uZp*+ imhl f L p-> + (3-12) 

Again the constants hi and <po are arbitrary, but to obtain a solution that asymptotes to the 
behaviour of section IBTTj we have to choose a particular value of hi for any fixed value of j^. 



3.3 Numerical solutions 

The expansions above are essential in the analysis of solutions to the system ( 13. lft - (13. 4ft . 
but their knowledge is not enough to guarantee the existence of an overall solution that 
connects the given UV and IR behavior. We have checked that such connecting solution 
indeed exists for the three kinds of IR behaviours. In this section, we plot some sample 
numerical computations. 

Figure [2] gives some examples of the numerical integration for the type I IR behaviour. 
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H[p] G[p] 




P 



P 



Figure 2: Type I solutions obtained numerically. We have fixed C = 0, N c = 1 and defined 
p = as the point where if (G) reaches a minimum for Nf < 2N C (Nf > 2N C ). We plot 
the solutions for different values of jf: dotted, dashed, dot-dashed and continuous lines 
correspond respectively to Nf = 1.2, 1.8, 2.2, 3 N c . 

As we have already anticipated, requiring that the type II and III expansions are connected 
with the UV asymptotics of section 13.11 imposes an additional condition, fixing the value of 
the integration constant hi in terms of C and jf. This additional condition produces the 
plots in Figure [Q Given that we choose the suitable value of hi, it is then possible to 
numerically integrate equations (13. ip - (13.41) . and find a complete solution for the type A 
backgrounds. For type II and type III solutions, we present examples of these solutions in 
Figures [3] and HI for different values of C and j^. 

3.4 UV expansions revisited 

The reader may wonder how it is possible that different IR behaviours asymptote to the 
same UV when it seems that the UV expansions (13.51) . (13.61) have no free parameters apart 
from C. Indeed, those expansions are not complete and there is an extra free parameter in 
the UV if one considers the possibility of exponentially suppressed terms: we can write, 

H = Q {p) + e- 4p p' a Qi(p) + 0(e- 8p ) (3.13) 
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H(p) H(p) 




Y(p) Y(p) 
I 1 7 . 5 1 




P 



P 



Figure 3: Type II solutions derived numerically. In these plots we present the results for C = 
— 1 and different values of ^f: dotted, dashed, dot-dashed and continuous lines correspond 
respectively to Nf = 1.2, 1.8, 2.2, 3 N c . The graphs on the left-hand side represent the 
solutions for the background on a large range of p, while the plots on the right-hand column 
focus on the behaviour of the functions close to the origin. Tilded function are rescaled by 
N c as H = jr, etc. 
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Figure 4: Type II and III solutions derived numerically. In these plots we present the 
behaviour of type III solutions, and compare it with small-C type II solutions. Clearly the 
transition through C = is almost continuous. Here we have fixed Nf = 1.5N C . Dotted 
lines correspond to the type III solution, while dot-dashed and continuous lines correspond 
to type II backgrounds with respectively C = —0.1 and C = 0.1. 
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where Q are the expansions (13.51) . (13.61) with: 

2N C -3N, 

3(2N C -N f ) y i ' 

a =w0h W>2JVJ (314) 

The Qi(p) are polynomials such that Qi(p) = Qi$ + Qi,ip _1 + Qi,2P~ 2 + Qi,o is a free 

parameter and the rest are fixed in terms of it. We have, thus (once C, N c , Nf are fixed), a one 
parameter family of UV solutions, labelled by Qi,o- Different choices Qi,o should modify the 
IR behaviour and, in particular, make the system fall into one of the qualitatively different 
IR possibilities. However, checking this numerically requires great numerical precision and 
we could not bear out this expectation explicitly. 

It would be interesting to know the physical interpretation of this parameter Qifl. Since it 
comes multiplying a term with e~ 4p , it must be related to the vacuum expectation value of an 
operator of UV dimension 6, here we anticipate the UV radius-energy relation log[/i/A] = |p 
that we will use in equation (14.61) . The simplest guess is an operator of the type {QQ) A (since 
in the UV 7q — > — |, at least when Nf < 2N C ). A single vev of this type would spontaneously 
break the R-symmetry, apparently in contradiction with the described solution. However, 
since there are infinitely many smeared flavored branes in the construction, it is natural to 
think that, if in fact such vevs are turned on, they may also be smeared effectively restoring 
the U(1)r. This is similar to the discussion in [13] where, considering smeared masses 
for hypermultiplets, the R-symmetry that would be broken for a single mass is effectively 
restored. 



3.5 A comment on integration constants 

The different expansions are written in terms of some integration constants which are related 
among them. In order to avoid confusion, we clarify here the number of parameters on which 
the different families of solutions depend. 

We have four first order equations (I2.2l) - (l2.5p . Let us fix Nf/N c (notice that, mantaining 
this ratio fixed, N c can be reabsorbed by simultaneously rescaling H, G, Y). Then a solution 
of the system depends on four integration constants. One is an additive constant for that 
does not affect the integration of the other functions. Since the BPS eqs do not depend 
explicitly on p, a second integration constant is just a shift of p, which we can fix without 
loss of generality. In type I, we can define p = as the point where H' (or G') vanishes, as 
in figure [2l For type II and III, we have defined p = as the point where H (and/or G) 
vanish. The third integration constant is fixed by requiring that the system goes to the linear 
dilaton UV (so there are not growing exponentials in H). From the IR point of view, this 
is, for instance, what fixes h\ in terms of C in the type II IR expansions, see the discussion 
in section 13.2.21 
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Thus, for fixed Nf/N c and fixing also the additive constant of 0, we are left with the fourth 
integration constant, which gives a one-parameter family of type I and type II solutions. For 
type III, where we further require C — 0, there is a single solution (a limiting point of the 
type II family). 

After these discussions on the solutions, we will focus on obtaining physics predictions 
from them. To this we turn now. 



4 Physics predictions of the new solutions 

In this section, we study some non-perturbative aspects that can be learned from the solu- 
tions presented in section [3J We will discuss expressions for the gauge coupling and anoma- 
lous dimensions of the quark superfields, give predictions for their values and compute the 
beta function in the UV of the theory (matching results known from field theory). The pre- 
dicted anomalous dimension for the quark superfield will be in agreement with the results 
of eq. ( 12.191) and those spelled out in section 12.3. 1L We will present similar results for the 
beta function of the quartic coupling (for which we will provide a geometrical definition). 
We will study the Wilson and 't Hooft loops and address in detail the issue of (continuous 
and discrete) anomaly matching. 

Let us remind the reader that we take units so that a' — g s — 1. With this choice, the 
tension of a D5 brane is T 5 = -^ys- 

4.1 Theta-angle and gauge coupling 

We will give an expression for the theta-angle and gauge coupling that can be simply derived 
by considering the action for a D5 brane on the submanifold (t,x,9 = 9, tp = 2n — (p), at 
constant radial coordinate. This is a well known definition, originally proposed in [35]. In 
[IT] we defined the gauge coupling using an instanton, whose action is computed using a Dl 
brane wrapping the two-cycle {9 — 9, (p — 2tt — (p). This two-manifold is special, being the 
only shrinking closed manifold in the geometry. As expected, the results using D5 or Dl 
probes do coincide. 

We consider the dynamics of a D5 with gauge fields excited on the R 1 ' 3 directions, 
obtaining (in string frame), 

S D5 = -T 5 J d 6 xe-^-det(g + 2irF) + ^(2tt) 2 T 5 J C 2 A F 2 A F 2 (4.1) 

The induced metric and two-form on the probe D5 brane read, in string frame, 



r/s 2 — 



ind 



dx\ 3 + (H + G) (d6 2 + sin 2 9d<p 2 ) 



(C 2 ) ind = 2Nc - Nf (ij-iJo) sin 9d9 A d^ (4.2) 
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Following the procedure of [36], we can expand (14.11) up to second order in the gauge fields 
(this is the QFT part of the a' expansion), perform the integrals over the two-cycle and 
compare to the usual Yang-Mills action in order to obtain the couplings. 
For the theta-angle, we get: 

e.(!^z> (,_,„,. (4 ,) 

Since changes in are associated with traslations ip — > ip + 2e then only some values of 
are allowed and they correspond to changes in —>■ + 2kn, which nicely 



2kr, 



2N C -N f 

matches what we wrote about anomalies between eq. (I2.23j) - (I2.25[) . 
The gauge coupling in our SQCD theory takes the form, 



8vr 2 



2(H + G). (4.4) 



For type I solutions, the series expansions of section 13.2.11 show that the gauge coupling 
vanishes in the IR. The expansions of section 13.2.21 for type II solutions show a version 
of 'soft confinement' (with bounded gauge coupling in the IR), while for type III more 
conventional confinement is obtained (with divergent gauge coupling in the IR). Figure [5] 
shows the behavior of the coupling as a function of the radial coordinate. 

4.1.1 Beta Function and Anomalous Dimension 

We can compute the beta function of the gauge coupling as defined in (14. 4ft . As we stated 
below eq.1 12.15p . we will be interested in the wilsonian beta function, that takes the expres- 
sion, 

-p- <91og(/i/A) 

Now, we need a relation between the radial coordinate p and the scale p in order to perform 
the derivative in (14.51) . In section 5.5 of [UJ we proposed that the radius-energy relation, in 
the UV of the gauge theory, that is for large values of the radial coordinate, is: 

log[/VA] = \p , (UV, p - oo) (4.6) 

This relation was proved in the unflavored backgrounds by identifying the gaugino condensate 
(an operator of dimension 3) with a deformation of the background that asymptotically 
behaves as e~ 2p , see [36, 37J. Since in the flavored backgrounds, the same kind of deformation 
(leading to the type N solutions) still has the same exponential behaviour [IT] . (14. 6p is 
trustable in the present setup. 
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Figure 5: These plots represent the behaviour of the gauge coupling as defined in equation 
f)4.4p in type II (top) and type III (bottom) backgrounds. For the top plots we have taken 
C — — 1 and different values of ^f. The correspondence between line styles and ¥f is the 
same as in figure [3) For the bottom plots, we have taken Nf = 1.5N C and C = —0.1 and 0.1 
(type II), and C = (type III). For the correspondence between values of C and line styles 
we refer the reader to the legend of figure HI The plots on the right are a zoom around p = 
of the plots on the left. 



We now proceed doing the computation (14. 5 p in the UV (for p — > oo). After substituting 
f )4.4p . (14.61) and (12. 2p . (12. 3p . one obtains an expression for the anomalous dimension of the 
quark superfield 7q in terms of the background functions, 

Inserting (13.51) . we find the UV behaviour of the anomalous dimension for Nf < 2N C : 

1 127V 

7o c —^ , (N f < 2N C ) (4.8) 

1Q 2 32 (2 N c -N f )p 2 y f c> K ' 

This expression is quite nice, because it matches the fact that the anomalous dimensions 
of the quark superfield approaches 7q = — | to have the R-charge R[Q] = 1/2, see the 
discussion in section 14.81 
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One may wonder what happens in theories with Nf > 2N C (which are Seiberg dual to 
those with N c < Nf < 2N C ). Introducing in ( 14. 7ft the expansion of Y for Nf > 2N C , eq. 
(P , we find: 

7q ~*-6^- (N f >2N c ) (4.9) 

1Q 2 Nf 8(N f - 2N c )p 2 w cJ K ' 

So 7q > — ~, as anticipated in section [2.2.21 
4.2 The quartic coupling 

We will now analyze similar issues for the quartic coupling in the superpotential of (12. lip . 
Let us start with the beta function. 

4.2.1 Beta function of the quartic coupling 

Let us go back to the field theory expression of the beta functions as we discussed them in 
eq. (l2.18p . We have, for the (dimensionless) quartic coupling: 

& = ^WAj = * + ^ (4 - 10) 

Using the radius-energy relation (14 .6p and the value of the anomalous dimension computed 
in (14.71) . we find: 

So, in view of (12. 3p . we can identify: 

\ogk = c+^-, (UV, p ^ oo) (4.12) 

where c is some constant. 

Let us analyze what happens with the quartic coupling when Nf < 2N C and Nf > 2N C . 
Whether it is relevant or irrelevant depends on the sign of 1 + 2jq = -J-(4y — N c ). 

When N f < 2N C : 

1 + 2 -<^- ^-N,f 2< ° (4 ' 13 > 
so the quartic coupling is relevant: starts at some constant in the UV and grows towards 
the IR. 

When Nf > 2N C , 

l + 2 lQ = ^-(N f -2N c )>0 (4.14) 

so the quartic is irrelevant and diverges in the UV. Notice the consistency of these results 
with the field theory discussion in 2!). .38] . 
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4.2.2 Geometric interpretation 



In principle, one could expect that the quartic coupling could be read as data of the geometry, 
i.e. as some quotient of the volumes of the different cycles of the manifold. In that sense, 
(14. 12[) is not very satisfactory. 

Without a definite understanding, we speculate in this section about a definition of the 
quartic coupling that makes use of the geometry of the internal manifold. Since (I4.12p is 
trustable in the UV, we require that our guess behaves similarly in the large p region. We 
start with the Nf < 2N C in which the coupling is relevant and bounded in the UV: 

k = e ° +f *> 1 + m^wf 1 + ' ' ' ' {N > < 2Ay (4 ' 15) 

where we have conveniently fixed c (an overall factor) to make the expression simpler. 

We will now propose a definition for the quartic coupling k in terms of the geometry. First, 
let us notice that, for fixed p, we have one well defined two-cycle, a well defined one-cycle 
and a couple of well defined three-cycles. The two-cycle is the one mentioned before and 
given by 

S 2 = [6 = 6, <p = 2n-0\, Vol[S 2 ]=Ane^(H + G). (4.16) 

The one-cycle, along the worldvolume of the flavor branes, is parametrized by the angle ip 
and its volume is 

S 1 = Vol[S 1 } = 47ie^ /2 W. (4.17) 
It is possible to define two different 3-cycles in the internal manifold 

S 3 = (6,(p,i/>), Vol[S 3 } = (47r) 2 e 3 / 2< Wr, S 3 = (0,&-0), Vol[S 3 ] = (Anfe^GW 

(4.18) 

We define the quartic coupling as a quotient of the volume of the two-cycle times the one- 
cycle divided by the volume of one of the three-cycles. This definition makes sense since the 
coupling we are interested in (relevant for matchings between string and gauge theory) must 
be adimensional: 

Vol^]xVo^ 

Vol[S 3 ] 1 ' 

Then: 

G N c , . 

K = 1+ H = 1+ 2(2N C -N f / + - < 4 - 20 > 

When studying the Seiberg dual magnetic theory, one should change S* 3 — > S 3 according to 
the prescription (12.211) . This amounts to the interchange H <-> G. 

When Nf > 2N C , the quartic coupling grows exponentially towards the UV and the 
definition (I4.19p should be taken carefully. This is perhaps not surprising since the field 
theory description with a quartic term breaks down and needs some UV completion. We 
relegate further discussion to Appendix [A] 
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4.3 A nice matching with the Nf = 2N C case 



In this section, we will show a connection between our solutions for Nf < 2N C with the one 
for Nf = 2N C and the one for Nf > 2N C , and analyse it in view of the properties of the 
quartic coupling discussed above. We will need to use that all of the solutions are of the 
form (12. ip . with the expansions in eqs. fj3.5[) - fl3.6[) and that for Nf = 2N C an explicit solution 
is given by [TT] (£ is a free parameter valued in < £ < 4), 

We now consider the solution with Nf < 2N C in the far UV (p — > oo), and change the number 
of flavors towards Nf — > 2N~ with a particular scaling limit 

lim (2N C -N f )p^ oo, (4.22) 

p^oo, Nf— >2N~ C 

so we have that the functions in eq. (l3.5p behave as 



L 2 J )p^oo, G « -p (4.23) 

According to (14.121) . the quartic coupling takes some finite value at this point. 

We see that we can match this with the Nf = 2N C solution (I4.2ip if we take £ = 0. After 
matching, we move freely in the parameter £ (since in the Nf = 2N C case, £ represents a 
marginal coupling [II]), from £ = towards £ = 4. At the point £ = 4 the functions take 
the values 

iV 

G^oo, H = <p = p (4.24) 

and the quartic coupling diverges. We can match (I4.24p with the expansion for Nf > 2N C 
in eq. (13.61) if we scale similarly as above, 

lim (Nf — 2N c )p — > oo. (4.25) 

The interpretation is the following: since when Nf > 2N C , the quartic is irrelevant, we 
can only connect those theories with the end of the Nf = 2N C marginal line in which this 
coupling is infinite. On the other hand, when Nf < 2N C , the quartic is relevant so these 
theories approach the conformal Nf = 2N C point by the opposite limit of the marginal line, 
where the quartic takes its smallest value. At Nf = 2N C , we can tune £ (accordingly changing 
the quartic coupling) to move from one end of the line to the other. The process continuously 
connects, in the space of theories, the Nf < 2N C case to the Nf > 2N C one. 

This "flow in the number of flavors" provides a nice quantitative check and suggests that 
our discussion about the quartic coupling in section 14.2.11 is sensible. It would be nice to 
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understand if similar matchings and flows may happen for different solutions with Nf = 2N C 

Let us now study the potential between static quark- ant iquark and monopole-antimonopole 
pairs, using Wilson and 't Hooft loops. 



4.4 Wilson Loop and QCD-string 

In order to get some intuition of what will happen to the quark-antiquark pair when sep- 
arated, we can analyze a quantity associated with the tension of the putative QCD-string 
that forms between them. The functional form of this tension (that can be thought of as 
the tension of an Fl-string when taken to the IR of the background) is given by, 



Tqq ~ y/gtt9xx\iR — e \IR 

1 QQ ! TypeII ^ e 

TQQ t TypeIII ~^ (4-26) 

Tqq gives a qualitative way of evaluating the tension in the IR dynamics. The expressions 
(I4.26P tell us that type I solutions will have vanishing QCD-string tension (indicating that 
quarks are free), type II and type III solutions will show finite tension, indicating confine- 
ment. This is in agreement with the behavior of the gauge coupling in eq. (14.41) . Of course, 
the presence of dynamical quarks will produce screening (string breaking) when the energy 
stored in the QCD-string (14.261) is of the order of the lightest meson in the theory. 

Let us now compute the Wilson loop following the standard prescription of [39] . Consider 
the action for an Fl string parametrized by 

t = r, x = a, p = p(a) (4.27) 

with x taking values in [0, L] (such that p(0) = p(L) = oo) and r in [0,oo). The induced 
metric on this string is, 

ds 2 ind = eH - dr 2 + (1 + AYp ,2 )da 2 ] (4.28) 
So, the Nambu-Goto action for this string is 

S = -!- / dae\ll + AYp' 2 (4.29) 
2ira' J v 

This is in the form to follow the usual treatment for Wilson loops in [ID] , and write expressions 
for the (renormalized) energy E and the separation L. There is a one-parameter family of 
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string profiles depending on the minimal value of p reached by the string, which we call po- 



E{p c 



1 



2vra' 



r°° rpo 

^ po) L + 2 2v / F(v / e^-e 2 ^°) - e*)dp - 2 / 2V¥e*dp 

J po Jo 



These integrals can be studied numerically for our solutions. We will present an example in 
figure at the end of the next section. 

4.5 t'Hooft loop 

It is interesting to perform a qualitative analysis similar to the one around eq. fl4.26p . In 
this case, the effective tension for the monopole-antimonopole pair is given by the effective 
tension of the magnetic string, represented by a D3-brane wrapping the two-cycle (I4.16P and 
extended along (x, t) stretched at the bottom of the geometry. See eq. (14.321) and below for 
details, 

T mm ~ e *{H + G)\ m 

T m fh,TypeI * ^ N j 
-L mm,TypeI I ' C ^ 

Tmm,,TypeIII ~^ (4-31) 

T mm gives a qualitative way of evaluating the tension of the monopole-antimonopole pair in 
the IR dynamics Intuitively, eq. (14.3ip indicates that type I and type II solutions confine 
the monopoles, while type III solutions exhibit free monopoles. 

Let us now study in detail the t' Hooft loop that can be thought of as the "Wilson loop for 
magnetic monopoles" (the rigorous definition is in terms of the commutation relations with 
the Wilson loop operator). It is a similar computation to the one for the Wilson loop in the 
previous section, as we explained above, for a extended wrapped D3 brane on the cycle £4, 

£ 4 =(t = T, xi = a, p(a), 6 = 6, (p = 2ir-<p). (4.32) 

This gives an induced metric that in string frame reads, 

ds] nd = \ - dt 2 + (1 + AYp' 2 )da 2 + (H + G)(d6 2 + sin 2 6d<p 2 )] (4.33) 

and the action (per unit time) for the 'effective string' obtained after integrating out the 
two-cycle is, 



S = AnT 3 J da^f+gip' 2 , f = e*(H + G), g = 2fW 



(4.34) 
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We observe that this expression is "Seiberg invariant" according to the definition (I2.20p . and 
analogous formulas to those in eqs. ( 14.301) can be used. They read 



L( P o) 



E = 4vrT, 



f(Po) 



p0 f yj Pip) - PiPo) 



dp 



f(po)L + 2 



gjp) ( 
Hp) ^ 



P(p)-P(s )-f(p))dp-2 g(p)dp 



(4.35) 



An example of numerical analysis for the expressions (I4.35P using the functions (14.341) is 
plotted in figure El 





2 . 5 



3.5 
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Figure 6: As a typical example, we plot E vs. L graphs for Nf = 1.2N C . The dotted line 
corresponds to the type I IR behaviour, the dashed line to type II (with C — — 1) and the 
solid line to type III. The graph in the left is for Wilson loops as discussed in section 14.41 
while the one in the right corresponds to the 't Hooft loops introduced in section 14.51 



Our interpretation of figure [6] is as follows: first, it is interesting to notice that the graphs 
for the backgrounds with different IR behaviours share some qualitative properties, suggest- 
ing a smooth interpolation among the different solutions. In all cases, there is a maximum 
length reached by the loop which we interpret as string breaking due to pair creation. % 
Moreover, in both graphs, there is a line (the type III one for the Wilson loop and the one 
corresponding to type I for the 't Hooft loop) which stretches more than the rest indicating 
linear confinement (before the string is eventually broken). This qualitatively suggests a 
change of roles of quarks and monopoles between type I and type III backgrounds. We will 
elaborate on this in the next section. Finally, let us point out that there are also minimal 
lengths in the graphs. Presumably, this is due to the growing linear dilaton, which signals 
the failure of the supergravity description we are using at asymptotically large p. With such 
an UV completion, there should be a minimal length (L min ) of the field theory that can be 

5 In [41] . it was argued in a similar situation that such kind of behaviour did not correspond to string 
breaking. 
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reliably probed by the present computation. Indeed, the converging growing lines on the 
left of the 't Hooft loop graph are an artifact of this UV behaviour (length and energy in 
equations (I4.35j) grow with p when p is taken large). Thus, we discard this part of such 
lines as unphysical. The fact that the type III line of the 't Hooft loop graph only follows 
this unphysical behaviour indicates that in the corresponding field theory, monopoles are 
free (or they are screened at a length shorter than L min ). 

We close this section with a technical comment about the obtention of the graphs in figure 
|6j In order to compare the loops probing the different solutions, one has to ensure that they 
asymptote to the same UV (once we fix the additive constant of integration 0o for one of the 
solutions, the 0o of the rest have to be fixed accordingly). It is also necessary to take into 
account that having C ^ 0, the UV behaviour of H is different from the C = case (see 
eq (l3.5j) ). Thus, in order to compare with other solutions, the definition of the coordinate p 
has to be shifted p — ► p — 2 N 2 -N f ( so ^ ne or igi n °f space is not at p = any more). 

4.6 Back to the 1970's 

Briefly, we want to make contact with the material in section 12.41 

The intuitive analysis for the string tensions in eq. fl4.26p and eq. (14.311) . shows that in the 
three types of IR solutions the behavior for the potential between the (non-dynamical) quark- 
antiquark and monopole-antimonopole pair is quite diverse. Indeed, the intuitive analysis 
presented around eqs. fl4.26p - fl4.3ip shows that for small separations compared to the mass of 
a pair, type I solutions confine monopoles and let quarks move freely. The opposite occurs 
in solutions of type III, while type II solutions present a situation in which quarks and 
monopoles experience a growing force when separated 0. 

Another important observation is the fact that one can smoothly interpolate between the 
solutions of types I, II and III by changing the values of the integration constants. Indeed, 
one can see that in the limit of vanishing constant C in type II solutions (see section 13. 2. 21) . 
the expansions change into those of type III in section 13.2.31 (this was repeatedly stressed 
when plotting the functions). The connection with type I solutions in section [3.2.11 must 
be thought in a more heuristic way, by noticing that if the constants C and hi in the type 
II solutions were absent the expansions could continue past p = to smaller values of the 
radial coordinate and the leading term in the expansion would go like the type I solutions 
(this connection could be made more precise in terms of the UV constants of integration 
C, Qi,o, see section [3T4l but this is hard to see in the numerical analysis). 

This 'heuristic matching' resembles the fact that in theories with fundamental fields there 
are no phase transitions in moving between a confining and a higgsed situation [33]. The 
different VEV's of the operator of dimension six mentioned in section 13.41 is what should 

6 We emphasize again that for long distances pair creation gives place to string breaking and the law 
changes. 
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take the theory from one vacua (a particular string background) to another. Together with 
the behavior for Wilson and 't Hooft loops, this leads us to propose that the solutions of 
type I, II and III behave like the Higgs, oblique confined and confined regimes of the SQCD- 
like theory in eq.( 12.11l) . In the oblique confined phase, some dyons condense and become 
free. This should be achieved in the present setup by the vanishing tension of a wrapped 
extended D3-brane with some F-strings disolved inside it, i.e. with the worldvolume electric 
field turned on. We leave this computation for the future. 

The qualitative picture described above is quite appealing. However, we should write here 
a word of caution. The tensions ( 14.261) and ( 14.311) depend on a combination of functions 
at the singularity; even when the singular behavior cancels from the expressions of ( 14.261) 
and (I4.3ip -a manifestation of the 'good' character of the singularity- these results should 
be taken with care. It could happen that the criterion we are using for 'good singularities' 
is not enough and that not all the discussed solutions should be accepted as physical. 

None of the described behaviours corresponds to an IR conformal fixed point. We remind 
the reader that SQCD with a quartic superpotential flows in the IR to Seiberg's conformal 
fixed point of SQCD [29] (at least in the ratio Nf/N c appropriate to fall in the conformal 
window). We believe that this does not show up in our set-up since, as explained in section 
12.2.11 the theory dual to the string solution differs from the usual SQCD+quartic near the 
IR, presumably due to a non-trivial VEV. 

4.7 A comment on the running couplings 

Coming back to the expressions for the Wilson and 't Hooft loops in eqs . ( 14 . 291) - ( 14 . 341) . we 
observe that they are invariant under the Seiberg duality transformation (12.201) . Even more, 
one can easily see that the gauge coupling defined in eq. ( 14.41) is also invariant. This may 
puzzle the reader, who may expect that while the electric theory is strongly coupled in the 
IR, the magnetic theory should be IR free and vice versa. This is not what we observe; our 
solutions predict that both theories are strongly coupled (or both have weak g\ M coupling) 
towards the IR. This situation is not totally original, it is indeed what happens inside the 
conformal window in Seiberg's SQCD. In contrast, in our case, we do not have a conformal 
window, but we will see that having two theories with similar IR behavior is not ruled out. 

Let us discuss this issue using the expressions for the wilsonian beta function in eq. (14.51) 
and the values for the anomalous dimensions we obtained in eqs.( l4.8l) - (14.9l) for the Nf < 2N C 
and Nf > 2N C cases respectively. First, we introduce the notation: the theory e is the one 
with gauge group SU(N C ) and Nf flavors. We will assume that in the e theory Nf < 2N C , 
and the superpotential is, 



We will denote by m theory the one obtained by Seiberg duality from the theory e. The 




(4.36) 
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c • 



m theory will have gauge group SU(n c ) with Nf flavors, of course the relation is n c = Nf—N ( 
This, together with the fact that in the e theory we set Nf < 2N C implies that Nf > 2n c in 
the theory m. The superpotential in the m theory is (we will ignore, since it plays no role, 
the double trace term), 

W m = k m MM + -qMq (4.37) 
H 

where we have used that the new fields in the Seiberg dual theory are the meson My (not 
charged under the SU(n c ) gauge field w a ), the quarks q and the antiquarks q. 

Now, let us compute the beta functions in each theory. For this, we will need to know the 
anomalous dimensions of the fields. Our solutions predict (in the UV) anomalous dimensions 
for the e quarks 7q and the m quarks r y q given by eqs. fl4.8p and (14.91) 
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^ = - y --- 3 -^2- 6 iV>- (8iV / -16n c )p 2 - (4 ' 38) 

Let us neglect the small 0(p~ 2 ) correction in the far UV (p — > oo region), keeping only the 
leading terms. Now, we compute the beta functions. We use the values for the anomalous 
dimensions in eq. (14.381) . obtaining that the wilsonian beta functions for the two theories are 
given by, 

(3^2 = 3N C - N f (l - 7q ) = %N C - N f ), 
(3^2 = 3n c - N f (l - lq ) = -(Nf - 2n c ), 

(4.39) 

This shows that either both the electric and magnetic theories confine or they both are IR 
free. The result in eq. (14.391) is precisely the one we obtain if we use the UV expansions for 
H(p),G(p) in section [37T1 together with eqs . (14 . 71) - (14 . 9 [) . We now move to the matching of 
anomalies. 



4.8 Global Anomaly matching 

In a field theory like SQCD or ours in eq. (12. 1 lj) there are two types of currents, global 

(denoted below by j g ) and local or gauge currents, denoted by J L . Hence, there are four 
possible 3-point correlators of currents one may compute 

<J L J L J L >, (4.40) 

< JgJgJL >, (4.41) 

< j g JlJl >, (4.42) 

< JgJgJg > ■ (4-43) 
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The first and second correlators (I4.40p - fl4.4ip must be zero, or the gauge symmetry is anoma- 
lous and the theory is inconsistent. 

The correlator with two local currents in eq. ( 14.421) if nonzero, indicates the presence of an 
anomaly in the global current < d^jg >7^ 0. While for non-anomalous symmetries in the 
sense of eq. fl4.42p . the correlator for three global currents in eq. (14.431) must match in the two 
descriptions of the theory. This is the 't Hooft anomaly matching condition. 

We will start by studying the anomaly for the R symmetry current, as expressed by the 
correlator in (14.42[) . and related to the results of sect ions \2 . 3 . 1 1 and 14 . 1 I This is what appears 
in the changes for the angle. A general expression for this anomaly is, 

if) -> e iel5 i), A0 = e n r T(r)R[ij)} (AAA) 

r 

where the ip denotes the massless fermions of the theory, n r is the number of particles in 
the r representation of the gauge group, T(r) is the weight of the representation (we remind 
the reader that T(adj) = 2N C , T(fund) = 1) and R[ip] is the R-charge of the corresponding 
fermion. We need to assign R-symmetry charges to the superfields. The gaugino has R[X] — 1 
and we want to determine the R-charges of the fundamental superfields. They cannot be 
read directly from the supergravity solution, but they can be extracted indirectly from the 
G-angle transformation. Let us start looking at the electric e-theory (with Nf < 2N C ) . From 
(I4.44p . we read (using R[4>q] = R[4>q], required by the symmetry): 

A9 e = e e (2N c + 2N f R[i/> Q ]) (4.45) 

The values of e e that leave G e invariant (modulo 27r) are written in (12.250 and immediately 
suggest R[iI>q) = — §. Thus, we are led the following assignation (for the fermions in the 
multiplet) under the global symmetry SU(Nf)y x x U(1)r, for the electric e-theory: 

i/> Q = (N f , 1,-1/2), = (#/,-!, -1/2), A = (1,0,1) (4.46) 

This implies that for the scalars of the quark multiplets R[Q] = R[Q] = \. Notice that this is 
consistent with the existence of the quartic superpotential preserving the R-symmetry at the 
classical level. Let us also discuss the relation of this R-charge to the anomalous dimension 
of the quark field. At a conformal point 

dimO = -R[0) (4.47) 

Using this relation and the UV value of the anomalous dimension of the squark (14. 8p . we 
find R[Q] = ^dim[Q] = |(1 + \^q) = |, in agreement with the discussion above. The 
relation (14.471) can be used here because our theory has an asymptotic UV fixed point (the 
beta-function for the gauge coupling becomes trivially zero since the gauge coupling itself 
goes to zero and the one for the quartic coupling also vanishes due to 7q = — ~, see f l4.10p ). 
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We now turn to the dual magnetic m-theory with n c = Nf — N c colors such that Nf > 2n c . 
We denote the magnetic quarks with lower case q. In this case, (14.441) yields: 

A6 m = e m (2n c + 2N f R[^ q ]) (4.48) 

Comparison to (I2.25P again suggests R[ip q ] = — |, which we believe is the charge relevant 
to the backgrounds discussed in this paper. However, there is another logical possibility 
which we find interesting and we will discuss it in section 14.8. 11 Therefore, we consider the 
following charges for the magnetic m-theory: 

ii M =(N,xN,, 0,0), A = (1,0,1) (4.49) 

This implies that for the scalars of the quark multiplets R[q] = R[q] = | and for the meson 
multiplet R[M] = 1, again allowing the superpotential to preserve the R-symmetry. However, 
given (14.91) . one can see that the relation (I4.47H is not satisfied. This should not come as a 
surprise, because the quartic coupling grows in the UV of the magnetic theory, see around 
eq. (14.14l) and the UV theory is far from a conformal point. 

Indeed, now that the assignations in eqs. (14. 461) . (14.491) have been given, the next step is to 
check whether 't Hooft matching in the sense of eq. ( 14.431) holds. Since the R-symmetry is 
anomalous, the continuous global symmetry is SU(Nf)v X U{1) B - The triangles that should 
be matched are, 

<SU(N f ) v >, < SU(N f ) 2 v U(l) B >, < >, <U(1) B TT> (4.50) 

where the last triangle representes the 'gravitational anomaly' with two insertions of T^. 
Let us quote the results in the e theory with the values of eq. (14.461) . we denote d ABC = 
Tr(T A [T B ,T c ]), 

< SU{N f ) 3 v >= d ABC N f N c (l 3 + (-1) 3 ) = 

< SU(N f ) v U(l) B >= tr(T A T B )2N f N c (l - 1) = 

< [/(1)| >=iV / iV c (l 3 + (-l) 3 ) = 

< U(1) B TT >= N f N c (l - 1) = (4.51) 
while for the m theory, with values (14 .49 p . we have 

< SU(Nf) v >= d ABC N f n c (l 3 + (-1) 3 ) = 

< SU{N f fyU{l) B >= tr(T A T B )2N f n c (— - — ) = 

< [7(1)| > =iV/ n c [(^) 3 + (-^) 3 ] = 

N N 

< U{1) B TT >= Nfn c (— -) = (4.52) 
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We see that there is anomaly matching between the electric and magnetic descriptions, 
satisfying 't Hooft criteria. Up to here, the computation is identical to that for the usual 
Af = 1 SQCD without superpotential. In [32], it was shown in that triangles involving 
discrete symmetries should also be matched. It is therefore a non-trivial consistency check to 
compute the values involving the remnant non-anomalous R-symmetry Z2N C -N f = Z^.-2n c - 
We now explicitly check the matching of discrete anomalies between the electric and mag- 
netic description. Following the paper [42] . we can compute the values for the triangles 
involving the R-symmetry. There is a technical subtlety explained in [32] ■ We need to 
have integer values for the different charges to match discrete anomalies, so they have to be 
rescaled. Thus, following the prescription in [32], we consider a^ = Zw e -2Nt R-symmetry. 
The baryonic and R-charges are 'rescaled' with respect to the values in the e-theory eq. (14.461) 
for the fermion in the multiplet and read: 

Nr - N - (Nt - N) 

^ Q = ( Nf ,^—^,-l), ^ = (N f ,- { f — C> ,-1), A =(1,0,2) (4.53) 

while for the dual theory m, rescaling in the same fashion the values in eq. fl4.49p . we have 

V„ =(#/,—,-!), V* =(#/,-—,-!) 
n n 

i/>m = (N f xN f , 0,0), A = (1,0,2) (4.54) 

The integer n is defined as the greatest common divisor of N c and Nf 

n = GCD(N C , Nf) (4.55) 

One has to study the following triangles: 

<Z k TT>, <Z k U(l) 2 B >, <Zl>, <SU(N f ) 2 Z k > (4.56) 

We will consider the difference between the triangles in the electric and the magnetic theory. 
They should match up to multiples of k = 4N C — 2 Nf [32] ■ We obtain the following results, 

< Z k TT > e = 2N c N f {-\) + 2(N 2 - 1), < Z k TT > m = 2n c N f (-l) + 2{n 2 c - 1), 
A(e - m) = (4.57) 

which is a multiple of k = 4N C — 2Nf as required by [42J . 
For the second triangle involving three 'discrete' currents, 

< Zl > e = 2N c N f {-lf + 8(iV c 2 - 1), < Z\ > m = 2n c N f {-lf + 8(n 2 c - 1), 

A(e -m) = 6N f (2N c - Nf). (4.58) 
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Again, the difference is a multiple of k. For the triangle involving two baryonic and a 
'discrete' current, 



< Z k U{\f B > e = 2N c N f (^ Nf - Nc ^j (-1), < Z k U(l)% > m = 2n c N f {-l) (^j , 

A(e-m) = 2^ J ±(N f -N c )(2N c -N f ), (4.59) 

we obtain the correct matching modulo k; finally for the triangle involving the non-abelian 
global symmetry 

< SU(N f ) 2 Z k > e = -2N c N f < SU(N f ) 2 Z k >= -2n c N f 

A(e - m) = -2N f (2N c - Nf), (4.60) 

We see that the assignation of charges appropriately rescaled (I4.53p . (14.541) satisfies the 
discrete anomaly matching up to multiples of k = (4iV c — 2Nf) as explained in [12] . This 
nice match indicates that the value for the R-charge suggested by the string background is 
consistent with the field theory computations using the charges in eqs. (14.461) .( j4~49l) . 



4.8.1 A different charge assignation 

We again consider the charges for the electric theory written in (14.461) . but give different 
R-charges for the magnetic theory. Coming back to ( 14.481) . a second possibility to match the 
supergravity result (12.251) is to take R[ip q } — \ — 2^-. 

Our second assignation of charges for the magnetic theory under the symmetries is, thus 
(like before, we give the charge of the fermions in the multiplet): 

/ (AT Nc 1 Hc \ I (AT Nc 1 Hc \ 

^ = ( ^'^'2-V' ^ = (iV/ '-^'2-V 

iPm = (N f x N f , 0, -2 + 4-f ), A = (1, 0, 1) (4.61) 

The nice thing about this second assignation is that it satisfies the relation in eq. (l4.47p . as 
can be checked by using (14. 9p (although, as explained above, there is no reason to expect 
(I4.47P to hold out of a conformal point). The downside is that this implies that the coupling 
k m in eq. (l4.37p must be charged, so the R-symmetry is completely broken at the classical 
level in the magnetic theory. This could only correspond to the dual of an electric theory 
in which the R-symmetry is spontaneously broken and, thus, it does not seem relevant to 
the type A backgrounds discussed in this paper. We believe that this second assignation 
could describe a situation like what happens for type N background in eq. (l2.6p or those 
solutions found in [11], which will have the same UV as our type A backgrounds. Since in 
this case there is no remnant R-symmetry, an anomaly matching computation as the one in 
(14561) - ( 14601 does not apply. 
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In summary, to match in field theory the R-symmetry anomalies to the ones computed by 
the string solutions, we proposed values for the charges of the multiplets under the global 
and discrete symmetries. Notice that the R-charges are checked by the string background in 
the sense that the anomaly is encoded in the background via the potential C(2) (see section 
12.3. If) . Both charge assignations we have given are consistent with this computation. On 
the other hand, the charges under SU(Nf) and U(1)b written above are suggested by field 
theory considerations but not from the actual string solution. However, it is remarkable 
that the string solution indirectly knows about them (at least in the case of assignation 
(14 .46 p . (14 .49 D ) throught the discrete anomaly matching computation displayed in eqs. (14.561) - 
(14.60}) . The results in this section together with the previous checks, make us very confident 
of the interpretations that we have put forward for the string background. 

5 Summary, conclusions and future directions 

Let us briefly summarize the different topics elaborated in this paper. We have found a set 
of solutions for backgrounds of type A in eq. (12.1 j) . proposed to be dual to a version of Af = 1 
SQCD (12. lip . With the new solutions, we predicted the behavior for different field theory 
observables and phenomena like Seiberg duality, gauge coupling and its running, Wilson 
and 't Hooft loops, anomalous dimensions of the quark superfields, quartic coupling, its 
running and anomalies (continuous and discrete). Based on those string theory predictions, 
we provided a large number of matchings between field theory expectations and string theory 
results. Interestingly, our solutions seem to realize the Fradkin-Shenker result for the non- 
existence of phase transitions when moving between higgsed and confining phases of a theory 
with fundamental fields [33] . 

This work leaves some open questions for future research. Indeed, it should be nice, among 
other things, to understand if the same 'phase structure' occurs in type N backgrounds in 
eq. (l2.6p and the behavior of the quartic coupling in that case. The picture of Seiberg duality 
for type N backgrounds could be understood on a similar way as we did in this paper. The 
case of massive flavors, either for type A or N backgrounds also presents interest because it 
gives a resolution of the IR (good) singularity and its many applications. These and other 
topics will be dealt with in |24j . 

It should also be interesting to apply our treatment for quartic couplings to cascading 
theories. Sorting out other dynamical aspects of the backgrounds presented here, like for ex- 
ample, the spectrum of glueballs, mesons and baryons [43J is of main interest. The Veneziano 
expansion we are adopting indicates that there will be interactions leading to decays between 
the first two types of excitations and suggests ways to study baryons as solitons in the BI-WZ 
action. 

It should also be interesting to extend all this story to the type II A case, using backgrounds 
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like the ones in [44] . Extensions to lower dimensional field theories may prove interesting for 
condensed matter problems. Finally, it should be nice to think about possible high energy 
phenomenology applications of this line of research. 
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APPENDIX 

A Beta function for K, m 

The computation of the beta function for the quartic coupling in the magnetic case seems 
to be tricky, so, here we will write the results we get and our interpretation. 

If we Seiberg dualize the coupling n e in eq. ( 14.201) we get the magnetic version K m = 1 + 
We compute the beta function from string theory using the radius-energy relation eq. (14.61) 
and we obtain 



1 ^ d\og[K m ] 3 d 



, a H\ 3 (N f -n c ) ,. 



K „r ™ Slog(f) 28p °V Gl ^(N l -2n c )p> 
Now, let us look at the QFT results. The superpotential after Seiberg duality is given by, 

W m = kMM + zqMq (A.2) 
The classical dimensions for fields and couplings are, 

[q] = m, [M\ = m, [k] = m, [z] = 1 (A.3) 

where we used that the "meson" M is a fundamental field in the magnetic theory and not a 
composite of electric quarks. The R-charges of fields and couplings, obtained from the main 
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part of the paper eqs. (14.6 II) 



arc, 



1 71 71 71 

%] = (2- 2 ^)> R M = (*Jf ~ 1). fl[A]=4(l-2^) (A.4) 

and the anomalous dimensions from the main part of the paper (for the meson it is computed 
assuming that the relation between dimension and R-charge is the one of a fixed point), 

12 1 n r 5 n r 3 (Nf — n r ) n r . . 



We can now compute 

& = *[-l+7«] = ft(12-^-6), ft = i[7, + ^M]=-^< w ^-"U ) 



On the other hand, we could have considered the superpotential after integration out of the 
massive meson in eq. (1A.2l) 

W m = 5qqqq (A.7) 

where the double trace term is present but has no relevance here. In this case the beta 
function from QFT has the expression, 

/% = J[l + aT j = (.-12^) + | c fcg_ (A,) 

It seems that the string theory computation in eq. tlA.ip is reproducing the beta function of 
a combination of field theory couplings in eq. (1A.6l) and eq. (1A.8I) . Indeed, the R-symmetry 
neutral combination that makes this work is K m = (Sk)^ 1 . 
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